Spatial modes in waveguided parametric down-conversion 
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The propagation of several spatial modes has a significant impact on the structure of the emission 
from parametric down- conversion in a nonlinear waveguide. This manifests itself not only in the 
spatial correlations of the photon pairs but also, due to new phase- matching conditions, in the output 
spectrum, radically altering the degree of entanglement within each pair. Here we investigate both 
theoretically and experimentally the results of higher-order spatial-mode propagation in nonlinear 
waveguides. We derive conditions for the creation of pairs in these modes and present observations of 
higher-order mode propagation in both the spatial and spectral domains. Furthermore, we observe 
correlations between the different degrees of freedom and finally we discuss strategies for mitigating 
any detrimental effects and optimizing pair production in the fundamental mode. 



I. INTRODUCTION 

The process of parametric down-conversion (PDC) is 
widely used in quantum optics and quantum informa- 
tion as a source of photon pairs. These photon pairs 
are generated inside a nonlinear crystal when a pump 
photon decays into signal and idler photons; typically 
the signal photon is used to herald the availability of an 
idler photon that can be used in subsequent experiments. 
The distribution of the pump photon energy between the 
daughter photons is determined by momentum conser- 
vation or, using the language common in optics, by the 
phase-matching conditions. 

After some initial discussion, experiments have demon- 
strated that for PDC in bulk crystals the orbital angular 
momentum (0AM) of the pump beam is conserved [l|,[2|- 
Hence, by controlling the spatial mode of the interact- 
ing beams a new degree of freedom becomes available 
for quantum information processing. It can be used to 
generate mode-entangled [H, Q and even hyperentangled 
bi-photon states for quantum metrology [j or quantum 



information-processing applications [El,!^]- 

In recent years attention has fallen on PDC in waveg- 
uiding structures 0, H, 0, The resulting implicit 

spatial-mode control has significant benefits for photon- 
pair production. Due to the tight confinement of the 
fields inside the waveguide (WG) the effective down- 
conversion rate into useful spatial modes increases sub- 
stantially. This, in turn, allows such sources to be 
pumped at greatly reduced power levels while still achiev- 
ing high photon fluxes relative to bulk crystal downcon- 
verters. In the ideal case only the fundamental spatial 
mode is present, and the waveguide output is naturally 
suited for efficient coupling to fibers yielding a robust, 
high-brightness photon-pair source. 

However, real- world waveguide sources usually deviate 
significantly from this ideal. Typical waveguides support 
the propagation of several spatial modes for the interact- 
ing pump, signal, and idler photons which can assume 
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any of the spatial modes guided at their respective wave- 
lengths In the standard case of a roughly rectan- 
gular waveguide, the radial symmetry of the system is 
broken and hence the 0AM conservation condition is re- 
duced to that of parity conservation between the spatial 
modes of the three fields. Therefore 0AM is no longer a 
useful parameter. Instead the fields must be decomposed 
into the waveguide transverse field mode solutions. 



In addition to the modification of the spatial character- 
istics introduced by the waveguide, the spectral proper- 
ties of the PDC process are also affected. The energy dis- 
tribution between the signal and idler photons, and hence 
their wavelengths, is governed by the phase-matching of 
the longitudinal components of the wave vector. Every 
spatial mode corresponds to a different transverse — and 
therefore also longitudinal — momentum, so for every 
set of three spatial modes of pump, signal, and idler 
we observe photon pairs with a specific distribution of 
wavelengths. The resulting output from the waveguide 
contains many of these distributions coherently super- 
imposed on one another. For various applications this 
can be a significant problem if a heralded single photon 
in a single well-defined spatio-temporal mode is needed. 
Likewise second-harmonic generation (SHG), the reverse 
process of PDC, is affected. 



In this paper, we present a study of the impact of 
discrete spatial mode structure on PDC in a nonlinear 
waveguide. We consider the connection between the spa- 
tial and spectral degrees of freedom as well as both the 
challenges and benefits that arise from using a waveg- 
uided PDC photon-pair source in quantum information 
processing. In Sees. II and HI we present the theoreti- 
cal groundwork for dealing with discrete spatial modes in 
PDC. In Sec. IV we investigate the spatio-spectral mode 
structure of our waveguide utilizing SHG. We present ex- 
perimental results on the spectral impact of spatial multi- 
mode PDC and theoretical predictions in Sec. V. Section 
VI discusses our findings with respect to spatio-spectral 
correlations. 
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FIG. 1: (Color online) Schematic picture of waveguided PDC 
in a periodically poled rectangular waveguide. 



II. THEORETICAL BACKGROUND 

We analyze the process of PDC in the interaction pic- 
ture. In order to derive the two-photon state of the down- 
converted fields we follow an approach used in 0, • 
Assuming a nonlinear dielectric crystal with length L, 
a constant waveguide cross section A, and propagation 
along the z direction (Fig. [1]), the down-converted bi- 
photon quantum state 1?/^) (within first-order perturba- 
tion theory) is then given by: 
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The susceptibility x*^^^ (^) indicates the second-order 
nonlinear ity, V is the volume of the waveguide, and t' 
denotes the interaction time. The three interacting fields 
= V->^->'^) identify the pump field and the down- 
converted signal and idler waves, respectively. Note that 
we neglected the vacuum state which is of no interest in 
the scope of this paper. 

The boundary conditions imposed on the fields propa- 
gating inside a waveguide define a discrete and finite set 
of allowed transverse field modes. This is in contrast to 
the infinite, continuous set of transverse modes possible 
in bulk crystal PDC. We derive the joint spatio-spectral 
wave function of the photon-pair by quantizing the fields 
of the signal and idler beams in the waveguide basis [l5[ . 
For a fixed frequency we denote the spatial waveguide 
modes in terms of a discrete set of orthonormal mode 
functions Un{r) {n = 1, . . . , A/'), where r = (x, y) is the 
position in the transverse direction and z is the propa- 
gation direction of the pump beam. The vector Qn = 
{kx^ky) indicates the transverse momentum associated 
with the nth spatial mode of the field, and the reduced 
longitudinal wave numbers of pump, signal, and idler 

beams are given by = \J [uj n{uj) / cf — \qn\^ - 

Their positive and negative frequency parts can now be 
written as: 

X exp [t [f3irH^^.)z - co^t] } a(f)t(c,^) 
where fi = {s^i) labels the signal and idler fields and the 



sum embraces all spatial modes propagating inside the 
waveguide. In Eq. ([2]) all constants and slowly varying 
field amplitudes have been merged into the overall pa- 
rameter B' ^ and the operator aii^'^\uj^) corresponds to 
the creation of one photon with a given frequency in one 
discrete spatial mode. 

If we treat the electric field of the strong pump 
beam classically with the approximation of an undepleted 
pump it can be written as: 



Ep{r, z,t)=y2 a(w)4''^w|^^ (r) 

I Jo 
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+ C.C.. 



In the general case of ultrashort pump pulses the laser 
source can exhibit a broad-band frequency structure 
which is modeled by the pump function a{u;). In the 
following the field amplitude of each pump spatial mode 



(r) is denoted by An (in the experiment this may 
be tuned by changing the coupling into the waveguide). 

Using Eqs.([l])-(|3]) and following Grice and Walmsley 
[l^ we find for the two-photon state: 
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with B defining an overall constant. The indices Imn 
label the triplet of the pump mode /, signal mode m, and 
idler mode n. 

In Eq. (|4]) we introduced the value A/3imn{^s^^i) = 

PI^\uJs + ^i) - Pm\uJs) - f3n\uJi) - Pqpm, which cor- 
responds to a momentum mismatch between the dif- 
ferent propagation constants. This is corrected by the 
quasi-phase-matching vector i3qpm arising from a pe- 
riodic variation in the x^^^ nonlinearity in the z direc- 
tion fabricated in the waveguide to achieve perfect phase- 
matching [1^] . The spectrum of the down-converted pho- 
ton pairs is given by the joint spectral amplitude (JSA) 

flmn{^s,^i) = a{iJs^U;i)X(t)imn{^s,^i), whcrC ^(cJs +CJi) 

is the pump distribution and (j)imn{'^s^^i) is the phase- 
matching function. 

Note that Eq. @ is similar to the two-photon state of 
collinear, plane- wave PDC, neglecting the spatial struc- 
ture of the propagating waves 0|. However, by ex- 
plicitly considering the spatial modes propagating inside 
the waveguiding material, we can see that the generated 
biphotonic state is emitted into a superposition of inter- 
acting mode triplets {Imn). Each triplet typically ex- 
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FIG. 2: Different spatial modes excited in the down- 
conversion process lead to a superposition of Gaussian shaped 
frequency distributions in the generated two-photon state, (a) 
pump function, (b) phase-matching function, and (c) JSA. 



hibits a different overall down-conversion generation ef- 
ficiency due to the overlap between the three interact- 
ing fields Aimn- Moreover the triplet possesses a unique 
spectrum fimn{^s-,^i) because of the different longitudi- 
nal wave vectors satisfying the phase-matching condition 
A/^/mn(^s, ^i) = 0. In the experiment the spread of the 
photons into different spatial modes can be tuned by con- 
trolling the preparation of the incoming pump wave into 
different waveguide modes . In conclusion the multi- 
mode PDC state is represented as: 



\^) =B A^j'^Aimn j duJs duOi fi 
Imn ^ 

xa(«)tK)a«t(^,)|0,0) 



(5) 



The JSA of a typical two-photon state with only three 
different mode triplets is plotted in Fig. [2l As can be 
seen the differences in the momentum mismatch lead to 
slightly shifted phase-matching functions and the pho- 
tons are generated in a superposition of these three fre- 
quency amplitudes. Hence the down-converted photon 
pairs are emitted into a composite state entangled in fre- 
quency and spatial mode (or transverse momentum). 



III. WAVEGUIDE MODEL 

To model numerically the properties of the gener- 
ated two-photon state for a given waveguide architec- 
ture and down-conversion process, we assume a rectan- 
gular waveguide with width height a difference 
in the refractive index between the waveguide and sub- 
strate materials of An and a poling period A. The waveg- 
uide dimensions and the difference in the index of re- 
fraction are deduced from a measurement of waveguide's 
numerical aperture. The extracted parameters suggest a 
rectangular waveguide with parameters H = 6 /im, W = 
4 /im, and An = 0.01. In addition to the refractive index 
difference An between the waveguide and the substrate, 
our model also takes into account the air boundary at 
the upper edge of the waveguide. The waveguides have 
an effective poling period of A = 7.59 /im. The estimates 
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FIG. 3: (Color online) Microscope image of the waveguide 
cross-section. The height and the width (H xW — 6x4 fxm^) 
given by the numerical aperture measurement are in good 
agreement with the image. Note the different nomenclature 
for the biaxial crystal basis and the coordinate system used 
for the calculations. 
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FIG. 4: (Color online) The first four z' polarized spatial field 
modes at 800 nm propagating in our waveguide. 



for the waveguide dimensions were verified with a micro- 
scope as shown in Fig. [3l 

In our type-II PDC process y' polarized pump photons 
near 400 nm are down converted into y' and polarized 
signal and idler photons around 800 nm. The spatial 
modes propagating in this waveguide architecture can 
be calculated with a semi-analytical dielectric waveguide 
model [13]. The solutions of the transverse wave equa- 
tions correspond to sine and cosine functions followed 
by an exponential decay into the surrounding material 
(see Fig. H]). The spatial modes and their respective 
transverse momenta can be conveniently labeled by the 
number of nodes in the transverse x and y directions: 
{x nodes, y nodes). Our waveguide architecture shows 
multi-mode behavior at the spectral range of interest and 
supports spatial modes from (0,0) up to (3,5) and (1,2) 
at 400 and 800 nm, respectively. For given field modes 
we can readily evaluate both the propagation constant 
^Cp,s,0 _ i^jYi^ji^ ^Yid the coupling constant Aimn for 
each mode triplet in Eq. (|4|) (cf. Table [I]). 

In the special case of a rectangular waveguide the par- 
ity between the interacting modes is conserved (/i + + 
rrii = 2N^ A/" G IN) similarly to the conversion of 0AM 
in bulk crystal PDC. All down-conversion processes vi- 
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TABLE I: Five selected coupling constants for different down- 
conversion processes possible in our WG architecture. 
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FIG. 5: Frequency distribution of the two-photon state for 
the 60 most important mode triplets. The PDC spectra are 
spread over a range of more than 200 nm. 



dating the parity conservation exhibit a vanishing cou- 
pling coefficient Aimn- However perfect parity conser- 
vation would only be given if signal, idler, and pump 
waves propagate in the same mode basis, which in gen- 
eral is not the case due to their different polarizations 
and wavelengths. Hence weak coupling exists between 
mode triplets that do not conserve parity. 

Altogether we expect for this configuration that the 
down-converted fields will be emitted in a superposition 
of 720 different mode triplets each having a distinct JSA 
and coupling efficiency. With a 1.14 nm pump band- 
width the frequency spectra are spread over 200 nm in 
wavelength as can be seen in Fig. O 



IV. STUDY OF SPATIAL STRUCTURE BY 
SECOND HARMONIC MEASUREMENTS 

SHG is a useful tool to analyze the spectral effects 
of spatial-mode propagation in our waveguides because, 
unlike PDC, all three fields are intense. Hence it is easy 
to measure both their spectral and spatial characteris- 
tics. In this way we can anticipate the spatio-spectral 
structure of the PDC emission from the waveguide by 
studying its SH response. The rich modal structure of 
the potassium titanyl phosphate (KTP) waveguides has 




FIG. 6: (Color online) Schematic of SHG. The two gray (red) 
lines represent the two pump fields and the three black curves 
are three different phase-matching functions for three mode 
triplets. The black (blue) line at 45° is the locus of points 
picked out by the pump function as the fundamental wave- 
length is scanned. SH signal can be generated at the inter- 
sections of this line and the three phase-matching functions. 



been studied in the past by Roelofs et al.^ who reported 
the existence of higher-order spatial modes by investigat- 
ing frequency doubling |18| . 

In SHG, the signal and idler modes at the fundamen- 
tal frequency are combined to form a second-harmonic 
mode at double the frequency. Thus the wavelengths of 
pump and output are reversed and in analogy to PDC 
this process is phase matched at the same frequencies: 
an SHG signal will only be present at wavelengths where 
pump and phase-matching functions overlap. Tuning the 
pump (fundamental) frequency moves the pump function 
along a +45° line in the joint frequency space as shown 
in Fig. [6l As this pump function crosses each individual 
phase-matching function a SH signal will be observed. 
If we use a pump with a sufficiently narrow bandwidth 
and align the pump beam such that we excite multiple 
waveguide modes, we can probe the position of the phase- 
matching functions for each spatial- mode triplet. 

In our first experiment we pumped a 2.1-mm-long, 
type-II PP-KTP (periodically poled) waveguide with a 
narrow-band Tiisapphire laser [0.6 nm full width at half 
maximum (FWHM), 2.6 ps autocorrelation length, and 
~1.5 MHz repetition rate] and measured the second- 
harmonic response with a spectrometer (Ocean Optics, 
Mikropack) while tuning the pump from 780 to 820 nm. 
Simultaneously, we recorded the spatial structure of the 
frequency doubled mode with a beam profiler charged- 
coupled device (CCD) camera (uEye, IDS Imaging). We 
observed a SH response at four different wavelengths (Ta- 
ble [III)- The corresponding spatial modes can be identi- 
fied as (0,2), (0,0), (0,1), and (0,0) as shown in Fig. [71 
and the SH spectra of the different spatial modes lie 2-5 
nm apart. These modes propagate in the waveguide and 
could therefore be exploited to pump down-conversion 
processes. 

With the help of the SH signal we verified the results of 
our model fSec.lIIip, specifically that it correctly predicts 
the separation between different SH peaks. In order to 
match exactly the absolute wavelength values from our 
model with the measured values, we introduced a small 
correction into the pump momentum. This was used 
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FIG. 7: (Color online) Observed second- harmonic modes. 
The modes (1) (0,2), (2) (0,0), (3) (0,1) and, (4) (0,0) were 
recorded at wavelengths of 393.9 nm, 397.9 nm, 400.0 nm and 
404.6 nm respectively. (Axis labels correspond to the coordi- 
nate system in Fig. [3]) 
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FIG. 8: (Color online) Experimental setup for measuring the 
spectral marginals of the signal and idler. BBO, /3-barium 
borate crystal for SHG; DM, dichroic mirror; CF, color glass 
filter; HWP, half- wave plate; WG, waveguide; PBS, polarizing 
beam splitter; MM, multi-mode fibre; SM, single-mode fibre; 
SG, spectrograph with CCD-Camera. 
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TABLE II: Spectral correspondence of different spatial modes 
in SHG. The second column states the mode triplets that yield 
the recorded wavelengths. The third and fourth columns show 
the measured and predicted positions of the SH modes. 



as a fitting parameter to correct a discrepancy between 
the Sellmeier equations [l^ calculated for a rectangular 
waveguide and those for the real waveguide used in the 
experiment. It should be noted that the value of this 
correction was only around 0.4% of the pump momen- 
tum and was the only free parameter in the model. The 
measured and predicted positions of the SH peaks were 
then in very good agreement as shown in Table [III 

V. MEASUREMENT OF 
SPATIAL-TO-SPECTRAL COUPLING IN PDC 

In our next experiment we investigated the spatial 
structure of the signal and idler photons from PDC more 
directly by measuring the spectrum of the single counts in 
both signal and idler arms. Measuring the single counts 
of only one arm discards the events in the other. How- 
ever the corresponding spectral peaks of the two mea- 
sured down-conversion distributions can be connected by 
energy conservation. This type of spectral marginal mea- 
surement corresponds to an integration over one spectral 
dimension of the JSA. 

We pumped our down-converter with the SH of a 



TABLE HI: Spectral correspondence of different spatial 
modes observed in the marginal measurements of the PDC. 



mode-locked Tiisapphire laser (798 nm, 10 nm FWHM, 
and 250 kHz repetition rate). A fairly narrow-band pump 
was required in order to resolve the spectral peaks aris- 
ing in the marginal distributions from different spatial 
modes: hence the SH had a spectral width of 1.1 nm 
centered around 399 nm. Our nonlinear medium was a 
3.5-mm-long type-II PP-KTP waveguide similar to the 
one studied in Sec. [iVl Signal and idler were separated 
at a polarizing beam splitter; one arm was coupled into 
a single-mode (SM) fiber, while the other was coupled 
into a multi-mode (MM) fiber. The fiber outputs were 
connected to a high-sensitivity spectrograph (Andor) as 
shown in Fig. [8l By rotating a half- wave plate before the 
polarizing beam splitter we could direct either signal or 
idler beam to the MM or SM fiber. 

Coupling to MM fibers allowed the spectra of all of 
the different PDC spatial modes to be observed simul- 
taneously. The large number of spectral peaks indicates 
that several spatial modes were excited inside the WG 
(Fig. [9|). As expected, the peaks in the multi-mode sig- 
nal and idler spectra could be accurately paired according 
to energy conservation. When coupling to SM fiber the 
main spectral peak corresponding to (0,0) mode in signal 
(790 nm) and idler (806 nm) dominates indicating that 
only the fundamental spatial mode is efficiently coupled. 

A crucial factor in influencing the structure of the pho- 
ton pairs is the coupling of the pump mode into the 
waveguide. We examined the overlap between the in- 
cident free-space Gaussian pump mode and the corre- 
sponding pump modes in the waveguide basis, = 

dr uf^ (r) E^fJ^^ (r). For a Gaussian pump beam with 
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FIG. 9: (Color online) The observed noise-subtracted signal 
and idler marginals in black (blue) are in very good accor- 
dance with the predicted spectra in gray (red). 



the correct beam waist properly aligned with the center 
of the waveguide almost all of the energy is deposited into 
the (0,0) waveguide mode and very little is coupled into 
higher-order waveguide modes. Nevertheless given that 
small pump misalignments were present we estimated a 
relative coupling into each higher order pump mode of 
10% of the energy coupled into the (0,0) mode. Taking 
into account this coupling of the beam into the waveg- 
uide in conjunction with the conversion efficiency of all 
the interacting mode triplets and their respective phase- 
matching functions, we calculated the generated JSA. 
Hence we were also able to predict the measured marginal 
spectra with very high accuracy. We observed that the 
output is dominated by four main down- conversion pro- 
cesses (Table IIIT|): the coupling to other possible down- 
conversion processes is suppressed such that they cannot 
be seen above the noise level in the measurement. The 
remaining small deviations from the predicted marginal 
measurements can be explained by the different coupling 
efficiencies of the various modes into the MM fibers and 
the differences between the real waveguide investigated 
and our rectangular waveguide model. Note that the pre- 
dicted spectra of the SHG and the PDC measurements 
are in very good agreement with each other. Our model 
requires only one fitting parameter which is the small 
constant offset in the effective Sellmeier equations for the 
waveguide. 



VI. DISCUSSION 

Our results from the previous sections demonstrate 
that higher-order spatial mode propagation has signifi- 
cant consequences for photon-pair generation in waveg- 
uide architectures. The majority of quantum optics ex- 
periments require single photons in a well-controlled sin- 
gle mode in any given degree of freedom [2O, [2l| . Conse- 
quently, in most single photon experiments, it is desirable 
to suppress all higher-order spatial-mode processes and 
to promote the emission of the photon pairs into the fun- 
damental mode [(0,0) ^ (0,0) + (0,0)]. 

In order to achieve this with waveguided PDC, we first 
turn our attention to the measures one can take to pro- 
mote the coupling of the pump beam into the fundamen- 
tal waveguide mode. The pump beam should be spatially 
filtered to have a Gaussian profile, either by the use of a 
pinhole and lens pair or a single- mode fiber. This mode 
should be carefully matched to the waveguide diameter 
either by choosing the appropriate lenses for efficient cou- 
pling into the waveguide or by directly butt-coupling a 
single- mode fiber tip to the waveguide. In general a spot 
size with a diameter larger than the waveguide profile 
provides pump propagation that is closer to single mode, 
whereas a smaller beam will tend to exhibit more multi- 
mode effects. Fine tuning the coupling of the pump beam 
into the waveguide has a large impact on the modal struc- 
ture of the pump and the optimum position can be found 
by monitoring the PDC spectrum. In our experiments we 
have found that this careful pump coupling can very effec- 
tively suppress all down-conversion processes originating 
from higher-order pump modes. 

Having implemented these changes, any spectral peaks 
observed in the output spectrum beyond those of the 
fundamental modes originate from the down-conversion 
of the pump photons from the fundamental mode into 
higher-order signal and idler modes, as already shown 
in Fig. [9l Due to the different spectra connected to 
the different spatial modes, the output can be filtered in 
the spatial or the spectral domain to ensure collection 
of emission from the fundamental signal and idler mode 
only. By coupling the output beams into single-mode 
fibers one can affect a spatial-mode filter with relatively 
high contrast. The impacts on the frequency spectrum of 
such a spatial mode filter can clearly be seen in Fig. [10] 
where the coupling of the down-converted photons was 
changed from MM to SM fiber. This demonstrates that 
the spectra of the beam can be manipulated by opera- 
tions in the spatial domain. On the other hand, a fre- 
quency filter may be used to operate on the spatial struc- 
ture of the photon pairs, and could be adjusted to only 
transmit photons with spectra originating from funda- 
mental mode propagation. 

Alternatively this problem can be addressed in the pro- 
duction process by reducing the size of waveguide cross 
section to only guide signal and idler beams in the fun- 
damental mode. However since the wavelength of the 
down-converted pairs is approximately twice that of the 
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FIG. 10: (Color online) Noise subtracted SM fiber measure- 
ment (gray, red) and MM fiber (black, blue) marginal mea- 
surements. In these graphs the impact of spatial filtering on 
the spectral measurements can clearly be seen. The SM fibers 
select the signal and idler photons in the ground mode and 
suppress the side peaks from higher order spatial mode signal 
and idler photons. 

pump, the light inevitably propagates not only in the 
fundamental waveguide mode but also in several higher- 
order spatial modes. Therefore multi-mode propagation 
effects still have to be taken into account, yet the number 
or the excited different down-conversion processes will be 
significantly reduced. 

It is conceivable that the high down-conversion effi- 
ciency possible with a waveguide source in addition to 
these additional degrees of freedom may be used to cre- 
ate a multiplexed photon-pair source by simple spatial 
or spectral filter operations on the generated biphotonic 
states. On the other hand these modes could be uti- 
lized as an information-processing degree of freedom and 
hence achieving control over the spatial emission profile 
is of paramount importance. For example, the spatial 
correlations may be applied in the generation of hyper- 



entangled states offering new possibilities in the produc- 
tion of Bell states and additional robustness in quantum 
error correction due to the enlarged Hilbert space. 



VII. CONCLUSION 

In this paper we have derived the two-photon state gen- 
erated in waveguided parametric down-conversion, ex- 
plicitly taking into account the spatial mode propagat- 
ing inside the waveguide material. We have shown how 
the propagation of a multitude of pump, signal, and 
idler modes affects the spatial and spectral structure of 
photon pairs emitted from a nonlinear waveguide. We 
have observed these effects both through measurements 
of the spectra of second-harmonic light generated in a 
KTP waveguide as well as in marginal measurements of 
down-converted photon pairs. The spectra derived from 
the multi-mode treatment of the SHG and the down- 
conversion fit the measured data extremely well, yielding 
the possibility to precisely engineer further experiments. 
Furthermore we have shown that we can influence the 
spectra of the generated biphotonic states by operations 
in the spatial domain. 

We suggested how these effects may be controlled to 
give a truly single-mode photon source in the spatial do- 
main by either careful experimental design or waveguide 
engineering. These findings are of significance for quan- 
tum information experiments where nonlinear waveg- 
uides are used to generate photon pairs and may be har- 
nessed to generate hyperentangled photon states or may 
be applied to design multiplexed photon-pair sources. 

Note added: Recently we became aware of an indepen- 
dent theoretical examination of spatial entanglement in 
waveguided PDC [22[. Furthermore a recent experimen- 
tal investigation of higher order spatial modes in KTP 
waveguides with sum frequency generation was brought 
to our attention [2^. 
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